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Abstract 



We design an adaptive finite element method to approximate the so- 
lutions of quasi-linear elliptic problems. The algorithm is based on a 
Kacanov iteration and a mesh adaptation step is performed after each 
linear solve. The method is thus inexact because we do not solve the 
discrete nonlinear problems exactly, but rather perform one iteration of a 
fixed point method (Kacanov), using the approximation of the previous 
mesh as an initial guess. The convergence of the method is proved for 
any reasonable marking strategy and starting from any initial mesh. We 
conclude with some numerical experiments that illustrate the theory. 

Keywords: nonlinear stationary conservation laws; adaptive finite ele- 
ment methods; convergence. 

1 Introduction 

In this paper we consider quasi-linear elliptic partial differential equations over 
a polygonal/polyhedral domain C R d (d = 2, 3) of the form 



where a : fix R + — > M + is a function whose properties will be stated in Section 2 
below, and / G L 2 (fl) is given. These equations describe stationary conserva- 
tion laws which frequently arise in problems of mathematical physics [15]. For 
example, in hydrodynamics and gas dynamics (subsonic and supersonic flows), 
electrostatics, magnctostatics, heat conduction, elasticity and plasticity (e.g., 
plastic torsion of bars), etc. Some of these examples are better modeled by vari- 
ational inequalities (see [6] and the references therein), but these fall beyond 
the scope of this article, which attempts to set a first step towards understand- 
ing the convergence and optimality of inexact Kacanov-type iterations, in the 
context of adaptive finite element methods. 

For a summary of convergence and optimality results of AFEM we refer the 
reader to the survey [10], and the references therein. We restrict ourselves to 
those references strictly related to our work. 

Inexact adaptive finite element methods have been considered for Stokes 
problem in [1, 7] using Uzawa's algorithm. Briefly, a Richardson iteration is 
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applied to the infinite-dimensional Schur complement operator and in each it- 
eration, the elliptic problem is solved up to a decreasing tolerance. In [1] linear 
convergence is proved and in [7] the optimality of the method in terms of degrees 
of freedom is proved, after adding some new refinement steps to the algorithm. 

In [8], a contraction property is proved for an adaptive algorithm based 
on Dorflcr's marking strategy [2] for nonlinear problems of the type (1), using 
Orlicz norms to cope with the very mild assumptions on the nonlinear term a. 

In this work we impose stronger assumptions on a, which guarantee the 
convergence of Kacanov's iteration [6]. More precisely, we assume that a(-, •) is 
decreasing with respect to its second variable (cf. (17)), and fulfills condition (3) 
below; which are the same assumptions stated in [6], where they consider the 
iteration on a fixed space, either finite- or infinite-dimensional. Our focus is the 
convergence analysis of the adaptive method that results from performing one 
mesh adaptation in each iteration of the non-linear solver. This turns out to be 
a very realistic and efficient method, based on the sole assumption that a linear 
system is exactly solved in each iteration step. 

This paper is organized as follows. In Section 2 we present the class of spe- 
cific nonlinear problems that we study, and some of its properties. In Section 3, 
we present the inexact adaptive Kacanov algorithm and in Section 4 we state 
and prove the main result of this article, namely the convergence of the dis- 
crete solutions produced by the algorithm to the exact solution of the nonlincar 
problem. Finally, in Section 5, we present some numerical experiments which 
illustrate the theory, and explore the practical boundaries of applicability of the 
algorithm. 

2 Setting 

Let f2 C M. d be a bounded polygonal (d = 2) or polyhedral (d = 3) domain with 
Lipschitz boundary. A weak formulation of (1) consists in finding u G Hq(Q.) 
such that 

a{u;u,v) = L(v), Vu G H%(Sl), (2) 

where 

a(w;u,v)= / a( ■ , |Vw| 2 )Vu • Vv, Vu>, u, v G (fi), 
Jn 

and 

l(v)= f fv, yveH^n). 

We require that a is locally Lipschitz in its first variable, uniformly with respect 
to its second variable (cf. (15) below). On the other hand, we assume that a is 
C 1 as a function of its second variable and there exist positive constants c a and 
C a such that 

c a < a(x,t 2 ) + 2t 2 D 2 a(x,t 2 ) < C a , Va; G ft, t > 0, (3) 

where Diet denotes the partial derivative of a with respect to its second variable. 
The last condition is related to the well-posedness of problem (2) as we will show 
below. Also, it is possible to prove that (3) implies that 

c a < a(x,t) < C a , Vxen,t>Q. (4) 



2 



These same assumptions are stated in [6], where several interesting applied 
problems are shown to satisfy them. 

Additionally, it is easy to check that the form a is linear and symmetric in 
its second and third variables. Also, from (4) it follows that a is bounded, 

\a(w;v,v)\<C a \\Vu\\n\\Vv\\ n , Vw,u,v e ftf(fi), (5) 

and coercive, 

c a \\Vu\\l<a(w;u,u), Vw,u£H%(ty. (6) 
If we define A : Hq(£1) —> H^ 1 ^) as the nonlinear operator given by 

{Au, v) := a(u; u, v), V u, v £ H (Q) , 

then problem (2) is equivalent to the equation 

Au = L, 

where L E is given. Assumption (3) implies that A is Lipschitz and 

strongly monotone (see [15]), i.e., there exist positive constants C a and ca such 
that 

\\Au-Av\\ H - H a) <C A \\V{u-v)\\ Q , Vu^e^ffi), (7) 

and 

(Au-Av,u-v)>CA\\V(u-v)\\f l , Vu,veH^(Q,). (8) 

As a consequence of (7) and (8), problem (2) has a unique solution and is 
stable [14, 15]. 

3 Adaptive algorithm 

In order to define discrete adaptive approximations to problem (2) we will con- 
sider triangulations of the domain f2. Let 7o be an initial conforming triangu- 
lation of fi, that is, a partition of fl into d-simplices such that if two elements 
intersect, they do so at a full vertex/edge/face of both elements. Let T denote 
the set of all conforming triangulations of f2 obtained from 7o by refinement us- 
ing the bisection procedures presented by Stevenson [13]. These coincide (after 
some re-labeling) with the so-called newest vertex bisection procedure in two 
dimensions and the bisection procedure of Kossaczky in three dimensions [11]. 
Due to the processes of refinement used, the family T is shape regular, i.e., 

diam(T) 

sup sup =: kt < oo, 

TeT TeT Pt 

where diam(T) is the diameter of T, and px is the radius of the largest ball 
contained in it. Throughout this article, we only consider meshes T that belong 
to the family T, so the shape regularity of all of them is bounded by the uniform 
constant kt which only depends on the initial triangulation 7o [11] • Also, the 
diameter of any element T G T is equivalent to the local mesh-size Ht '■= IT] 1 /' 1 , 
which in turn defines the global mesh-size Ht '■= m&xHr- 
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For the discretization we consider the Lagrange finite element spaces con- 
sisting of continuous functions vanishing on dQ which are piecewise linear over 
a mesh T G T, i.e., 

V r := {v G H£(fl) | v w ^Vi{T), VTgT}. (9) 

We are now in position to state the adaptive loop to approximate the solution 
u of the problem (2). 



Adaptive Algorithm. Let 7o be an initial conforming 
of f2 and uq G be an initial approximation. Let 71 
and k = 1 . 


mesh 

= To 


1. u k :=SOLVE(u fe _i,T fe ). 




2. {r, k {T)} TeTk '■= ESTIMATE(ufc_i, Uk, %) ■ 




3. M k := MARK({i lk (T)} TeTkl T k ). 




4. r fe+ i := REFINE(7fc, Alfc, n) . 




5. Increment fc and go back to step 1. 





We now explain each module of the last algorithm in detail. 



The module SOLVE. Given the conforming triangulation T k of ft, and the 
solution from the previous iteration, the module SOLVE outputs the 
solution u k G Vfe := Y-f k of the linear problem 

a(u k -i;u k ,v k ) = L(v k ), V v k e Y k . (10) 

Remark 3.1 (Linear vs. nonlinear). Notice that while SOLVE requires the 
solution of a linear system, the usual discretization of (2) in consists in 
finding u k £ V k such that 

a(u k ;u k ,v k ) = L(v k ), V v k G Y k , (11) 

which is nonlinear. We propose here to make only one step of a fixed point iter- 
ative method to solve the nonlinear problem, and proceed to the mesh adapta- 
tion, whereas the usual approach would be to approximate the discrete nonlinear 
problem up to a very fine accuracy (pretending to have it exactly solved) before 
proceeding to mesh adaptation [5]. Each iteration of the adaptive loop is thus 
much cheaper in our approach. Our theory guarantees convergence of this algo- 
rithm, and the numerical experiments of Section 5 show that the convergence 
is quasi-optimal, although the latter is not yet rigorously proved. 

Remark 3.2 (Kacanov's Method). Let us consider problem (2) (resp. prob- 
lem (11) with k G N fixed). We denote the space Hq(SY) (resp. Vjt) by V, and 
the solution u (resp. u k ) by U . Given an initial approximation Uq G V of the 
solution U, we consider the sequence {Ui}i e ^ where U% G V is the solution of 
the linear problem 

aiUi-^U^v) =L(v), VdgV, i G N. 
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This is known as Kacanov's Method, and it follows [6] that the sequence 
{Ui}i£N converges to the solution U, provided D2d{x,t) < for all x G f2 and 
t > 0. 

Notice that our algorithm consists in performing only one step of Kacanov 
iteration in each step of the adaptive loop. 

The module ESTIMATE. Given T k , u k _\ and the corresponding output u k 
of SOLVE, the module ESTIMATE computes and outputs the local error 
estimators {r] k (T)}TeT k given by 

r, 2 k {T)-.= H 2 T \\R k \\ 2 T + H T \\J k \\ 2 dT , (12) 

for all T G 7fe. Here R k denotes the element residual given by 

R kw := -V- [a{- ,\Vu k - X \ 2 )Vu k ] - /, VTeT fe , 

and Jfc denotes the jump residual given by 

J k\ s : = 2 [M->|Vw fc _i| 2 )Vu fc )| r •n+(a(-,|Viife_i| 2 )Vu fe )| T , ■ n'] , 

for each interior side 5, and J k i s ■= 0, if 5 is a side lying on the boundary 
of f2. Here, T and T" denote the elements of T k sharing S, and n, n' are 
the outward unit normals of T, T' on 5, respectively. 

The squared sum of the a posteriori error estimators is an upper bound for 
the residual R(wfe) G 1 (SI) of which is dchned as 

(R(ufe),u> := a(u k -i;u k ,v) - L(v) = I (a( • , |V-« fe „i| 2 )Vit fc • Vji - /w) , 

Jo 

for u G iJg(ri). In fact, integrating by parts on each T G T k we have that 

(R(u fc ),«)= V ( f R k v+ f J k v), (13) 

TeTk \Jt JdT / 

and since (R(itfc),w/ £ ) = for v k G Vfe, using (13) and interpolation estimates, 
we arrive at the following upper bound: 1 

\{R(u k ),v)\ < Vk(T)\\Vv\\ Ukm , VvcH^n), (14) 
TeT fc 

where ui k (T) is the union of T and its neighbors in T k . 

The next result is some kind of stability result for the estimators, which is a 
property somewhat weaker than the usual efficiency, but sufficient to guarantee 
convergence of our adaptive algorithm (see also [12, 4]). In order to prove it, 
we assume that a is locally Lipschitz in its first variable, uniformly with respect 
to its second variable, as we mentioned before. More precisely, we assume that 
there exists a constant C a > such that 

\a(x,t) - a(y,t)\ <C a \x-y\, Vz, y G T, t > 0, (15) 

for each T G To- 

1 From now on, we will write a < b to indicate that a < Cb with C > a constant depending 
on the data of the problem and possibly on shape regularity of the meshes. 
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Proposition 3.3 (Stability of the local error estimators). Let {wfcjfcgN be the 
sequence of discrete solutions computed with the Adaptive Algorithm. Then, the 
local error estimators given by (12) are stable. More precisely, there holds 

Vk(T) < HVtifclLcD + II/IIt, VT £ Tk, 

for all k £ N. 

Proof. Let {u k } k <£N be the sequence computed with the Adaptive Algorithm. 
Let k £ N and T £ T k be fixed. On the one hand, using that Uk\T is linear, and 
thus Awfc = inside T, we have that 

\\R k \\ T = ||-V ■ [a(-, |Vti fc _i| 2 )Vu fc ] - f\\ T 

< ||V[q;(-, |Vw fc _!| 2 )] - Vw fc || r + . 

Since a is locally Lipschitz in its first variable (cf. (15)), it follows that if £ := 
Vwfc-i (constant over T), 



5j:<»fe|{| 2 ) 



<C a , VieT, l<i<d, 



and thus, 

ll^fc|lr^l|Vufc|| r + ||/|| T . 

On the other hand, if S is a side of T k shared by the elements T,T'eTt, 

HJfclls < ||K->|Vw fe -i| 2 )Vu fe ) k || s + ||(a(-,|V Ufc _i| 2 )Vu fe )| T ,|| s 
< 



|Vufc| T || s + ||Vu fc | T , | 
< ff- 1/2 ||Vu fe || T + i/ T , 1/2 ||V Ufc || T ' < H T 1/2 \\Vu k h 
where we have used (4) and a scaled trace theorem. Therefore, 



HT /2 \\Jk\\ dT <\\Vu k \\ UkiT) , 

which completes the proof. □ 

The module MARK. Based on the local error indicators, the module MARK 
selects a subset M. k oiTk, using any of the so-called reasonable marking 
strategies, such as the maximum strategy, the equidistribution strategy, or 
Dorfler's strategy [11]. More precisely, we only require that the set of 
marked elements Ai k has at least one element of Tk holding the largest 
local estimator. That is, there exists an element T™ ax £ Ai k such that 

Vk(T^) = maxr/ fc (T). (16) 

i t Ik 

This is what practitioners usually do in order to maximize the error re- 
duction with a minimum computational effort. 

The module REFINE. Finally, the module REFINE takes the mesh T k and 
the subset M k C T k as inputs. By using the bisection rule described by 
Stevenson in [13], this module refines (bisects) each element in Ai k at 



G 



least n times (where n > 1 is fixed), in order to obtain a new conforming 
triangulation Tk+i of O, which is a refinement of Tk and the output of 
this module. By definition, Tk G T for all k G N and the family of meshes 
obtained by the Adaptive Algorithm is shape-regular. Finally, it is worth 
observing that the resulting spaces are nested, i.e., Vfc C Vjt+i; this fact 
will be used in the proof of Lemma 4.2 below. 

Remark 3.4 (The adaptive sequence {itfc}fceN is bounded). By the coercivity 
of a given by (6), and using the definition (10) of Uk we have that 

||Vttfc||n < — a(uk-x; Uk,Uk) = —L(u k ) < - ^ H — — ||Vufc||n, 

Ca C a C a 

for all k G N, and then 

,, H^IU-i(n) 

Ca 

Therefore, {ufe}fcgN i s bounded in Hq(Q). 



4 Convergence Analysis 

In this section we prove the convergence of the sequence computed with the 
Adaptive Algorithm described in the previous section. We combine the ideas 
of the proof of convergence of adaptive algorithms for linear problems given 
in [9] and [12] with new techniques needed to overcome the difficulties arisen by 
the nonlinear nature of the problem treated in this paper, adapting some ideas 
from [3, 4]. 

As a first step to prove the convergence of the {wfe}fegN to the exact solution 
u, we prove that ||V(itfc — Mfe+i)||n — > 0, as k tends to infinity, for which we need 
the following auxiliary lemma. 

From now on we assume that 

a(x,ti) > a(x,t2) whenever < t\ < t<i and ie!!. (17) 

Lemma 4.1. Let us assume that a(x, ■) is a monotone decreasing function for 
all x G £1, i.e., (17) holds. Then 

J(v) - J{w) < -(a(w;v,v) - a(w;w,w)), Vv, w G i?o(0), 

where J{v) = J Q a(sv; sv, v) ds. 

Proof. Let v, w G Hq(Q). Then, a change of variables in the integral defining 
J(-) yields 

Iff /-IVi-l 2 HVuf \ 

J{v) - J(w) = - / I / a (x, t) dt- a(x, t) dt J dx 



nlV^I 2 

a{x, t) dt dx. 
in J\Vw\ 2 
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Since a(x, •) is a decreasing function for all x G Q, 
1 f f\ Vv \ 2 1 f 

- / / a(x, t) dtdx<- / a(x, | Vw| 2 ) (|Vif - |Vw| 2 ) da; 

= i (a(w; u, u) - a(w; w, w)) , 

and the assertion of the lemma follows. □ 

Lemma 4.2. Let {uk}keN denote the sequence of discrete solutions computed 
with the Adaptive Algorithm. Then, 

Urn \\V(u k - u k+1 )\\ n = 0. 

k—yoo 

Proof. Let {it/c}fceN be the sequence obtained with the Adaptive Algorithm. 
Since a is coercive (cf. (6)), and linear and symmetric in its second and third 
variables, we have that 

c a ||V(Mfe - u k+ i)\\n < a(u k ;u k - u k+1 ,u k - u k+1 ) 

= a(u k ;u k ,u k ) - 2a(u k ;u k+ i,u k ) + a{u k ; u k +i, u k +i). (18) 

From (10), since u k G Vfe+i, it follows that a(u k ;u k+ i,u k — Ufc+i) = L(u k —u k+ i) 
and thus 

a(u k ;u k+ i,u k ) = L(u k - u k+ \) + a(u k ;u k+ i,u k+ i). 
Replacing this equality in (18) and taking into account Lemma 4.1 we obtain 

c a \\y{u k - u k+ i)\\ 2 n < a(u k ;u k ,u k ) - 2L(u k - u k+1 ) - a(u k ;u k+1 ,u k+1 ) 

< 2J(u k ) - 2J(u k+1 ) - 2L(u k ) + 2L(u k+1 ) (19) 
= 2{F{u k )-F{u k+1 )), 

where J- := J — L, and therefore, {J-(u k )} ke jq is a monotone decreasing se- 
quence. 

On the other hand, {J-(u k )} ke pj is bounded below since 

F(u k )= I sa{su k ;u kl u k ) ds - L(u k ) 
Jo 

> 2 C all VM fc|| 2 2 - ||£|U-i(n)||Vu fe ||n > ^ — — • 

Finally, from the last two assertions it follows that {J-(u k )} k £f$ Q is conver- 
gent. Considering (19), we conclude the proof of this lemma. □ 

We show now that the sequence obtained with the Adaptive Algorithm is 
convergent, and more precisely, that it converges to a function in the limiting 



-irj(n) 



space Voo := l)Y k . Note that Voo is a Hilbert space with the inner product 
inherited from Hq(Q). 

Theorem 4.3 (The adaptive sequence is convergent). Let {u k } k ^ a be the se- 
quence obtained with the Adaptive Algorithm. Let 6 Voo be the only solution 
to 

Woo G Voo : a(Moo; Woo, v) = L(v), V II £ Voo- (20) 

Then 

u k — > Moo in Hq(Q,). 
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Remark 4.4. Notice that (20) always has a solution because (7) and (8) hold 
on Voo, which is itself a Hilbert space. 

Proof. Let {uk}keN be the sequence obtained with the Adaptive Algorithm. 
Let Moo € Voo denote the solution of (20) and Vk+i '■ Hq(Q) — > ^k+i be 
the orthogonal projection onto Vk+i- Since A is strongly monotone (cf. (8)), 
using (20) and (10) we have that 

C A \\V(u k - Uoo)||n < (Au k - AUoo^Uk - Uoo) 

= (Au k , U k - Uoo) - L(u k - Uoo) 

= {Au kl U k - Vk+lUoo) + (Au k ,V k +lUoo - Uoo) 

- L(u k - V k+ \Uoo) - L(V k+ iUoo - Uoo) 

= a(u k ]U k - U k+1 ,U k - "Pfc+lUoo) + (Au k ,Vk+lUoo ~ "oo) 

- L(Vk+lUoo - Uoo) 

< C a \\V(u k - Ufc+OHndlVufclln + ||Vuoo||n) 

)l|n, 

for all fc S No, where in the last inequality we have used (5). From Remark 3.4 
it follows that {ufc}fe e N is bounded in Hq(CI), and using Lemma 4.2, together 
with the fact that the spaces {Vfc}fc S N are nested and U ke n ^k is dense in Voo, 
we conclude that u k Uoo in Hq(£1). □ 

In order to show that the limiting function Uoo is, in fact, the solution of the 
problem (2) and thereby conclude that the adaptive sequence converges to the 
solution of this problem, we establish first two auxiliary results (see Lemma 4.6 
and Theorem 4.7). We need the following 

Definition 4.5. Given any sequence of meshes {Tk} k ^n C T, with T k +i a 
refinement of T k , for each fc € No, wc define 

T+ :={TeT k \TeT m , Vm>fc}, 7^:=Tk\T k + , 

and 

0+ := |J Uk {T), n° k := (J Wfc (T). 

T6T+ TeT° 

In words, T k is the subset of the elements of T k which are never refined in the 
adaptive process, and Tj? consists of the elements which are eventually refined. 

It can be proved [9] that if xq° denotes the characteristic function of 0^, 
then 

lim ||/ifcXn°IU«(n) = 0, (21) 

where h k £ L°°(f2) denotes the piecewise constant mesh-size function satisfying 
h k \ T := H T , for all T £ %. 

Since the error estimators are stable (cf. Proposition 3.3), using the conver- 
gence proved in the last theorem and (21) we can establish the following 

Lemma 4.6 (Estimator on marked elements). Let {{rjk(T)}T^T k } k<=N ^ e ^ e 
sequence of local error estimators computed with the Adaptive Algorithm, and 
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let {Mk}k<£N a be the sequence of subsets of marked elements over each mesh. 
Then, 

lim max r>k(T) = 0. 

Proof. Let {{Vk(T)}TeT k }keN a an< ^ {-Mk}keN be as in the assumptions. For 
each k £ No, we select Tk £ A4k such that rjk(Tk) — uiaxTeM k Vk(T). Using 
Proposition 3.3 we have that 

m{Tk) < \\Vu k \\ Uk{Tk) + \\f\\ Tk < llVufc-Vwoolln + HVuooll^TO + ll/lk, (22) 

where £ Voo is the solution of (20). On the one hand, the first term in the 
right hand side of (22) tends to zero due to Theorem 4.3. On the other hand, 
since T k £ M k CJi", from (21) it follows that 

\T k \ < |w fc (T fc )| < U d Tk < ||/ifcXno||i-(n) °> as fc ^ oo, 
and the last two terms in the right hand side of (22) also tend to zero. □ 

Using the upper bound (14), the stability of the estimators (Proposition 3.3), 
the facts that {ufe}fc S N is bounded (cf. Remark 3.4) and the marking strategy 
is reasonable (cf. (16)), and Lemma 4.6, we now prove the following important 
result. 

Theorem 4.7 (Weak convergence of the residual). If {ufe}fegN denotes the 
sequence of discrete solutions computed with the Adaptive Algorithm, then 

lim (R(ufc), v) = 0, for all v £ H£(S}). 

k-^oo 

Proof. We prove first the result for v £ H 2 (£l) n Hq(Q), and then extend it 
to H^fl) by density. Let p £ N and k > p. By Definition 4.5 we have that 
T p + C C Tk- Let Vk £ Y k be the Lagrange's interpolant of v. Since 
(R(v,k),Vk) = 0, using (14), and Cauchy-Schwartz's inequality we have that 

\(R(u k ),v}\ = \{R(v,k),v- v k )\ < £ Vk(T)\\V(v-v k )\\ Uk(T) 

rer k 

= E Vk(T)\\V(v-v k )\U {T) + J2 Vk(T)\\W(v-v k )\\ Uh(T) 
TeTj" TeT k \r^ 

< m{T p + )\\V{v v k )\\ nt + Vk (Tk \ T P +)||V(« - v k )\\ n0p , 

1 

where, for any 3 C Tk we hereafter denote (X^t<=~ 7 ?fc(^)) 2 by %(S). Taking 
into account Proposition 3.3 and the boundedness of the discrete solutions (cf. 
Remark 3.4) we have that r]k(Tk \ T p + ) < %(71-) < 1, and therefore, 

|(R(u fc ).u)| < (r) k (T p + ) + \\h p xa%h°°{(i)) K'Ih 2 (o), 

due to interpolation estimates. 

In order to prove that (R(uk),v) — > as k — >• oo wc now let e > be 
arbitrary. Due to (21), there exists p £ N such that 



\\hpXQ'>\\L^(Q) < 



s. 
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On the other hand, since C T k + C T k and the marking strategy is reasonable 
(cf. (16)), 

Vk(r p + ) < (#T/) 1/2 max % (T) < (#T/) 1/2 max m (T). 

Now, by Lemma 4.6, we can select K > p such that r/ k (Tjf) < e, for all k > K. 
Summarizing, we have proved that 

lim (RK), v) = 0, for all v € H 2 {VL) n H^(Q). 

k— >-oo 

Finally, since H 2 (il) n ^(17) is dense in TJq(JI), this limit is also zero for all 

veH^n). □ 

As a consequence of Theorem 4.7 we now prove that Uoo is the solution of 
problem (2). 

Theorem 4.8 (The limiting function is the solution). Ifu^ denotes the solution 
of (20), then is the solution of problem (2), i.e., 

a^oo^oo,") = L(v), Vji£H'(1]). 

Proof. Let ttoo be the solution of (20). If v € Hq(SY), and {itfc}fceN denotes the 
sequence of discrete solutions computed with the Adaptive Algorithm, then 

K«9o;woo,f) - L(v)\ = Kuo^Woo, v) - L(v) - a{u k ;u k+ll v) + a{u k ; u k +i, v)\ 

< (a^jUoct)) - a(u k ;u k+1 ,v)\ + \{R(u k+1 ),v)\ 

< ^(uoojMocf) - a(u k ;u kl v)\ 

+ \a(u k ;u k - Uk+i,v)\ + |(R(it fe+ i), v)\ 

< \\Auzo - Au k \\ H -i( n )\\Vv\\n 

+ C a \\V(u k - u fc+ i)||n||V«|| n + |(R(«fc),«)| 
<C A ||V( Uoo - Ufc )llo||V V ||n 
+ C , ||V(ujfc-ttfc + i)||n||V«|| n + KR(«fc),«)| 1 

where we have used that A is Lipschitz (cf. (7)) and a is bounded (cf. (5)). 
Using Theorem 4.3, Lemma 4.2 and Theorem 4.7 it follows that the right-hand 
side in the last inequality tends to zero as k tends to infinity. □ 

As an immediate consequence of Theorems 4.3 and 4.8 we finally obtain the 
main result of this article. 

Theorem 4.9 (Main result). Let {u k } k £jq denote the sequence of discrete solu- 
tions computed with the Adaptive Algorithm. If a(-, •) satisfies assumptions (3) 
and (17), then {u k } ke n converges to the solution u of problem (2). 

We conclude this section with a couple of remarks. 

Remark 4.10. The problem given by (1) is a particular case of the more general 
problem: 

f -V ■ [a( ■ , \Vu\ 2 A )AVu] = f infi 

u = on 90, 
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where a : ft x R + — > R + and / E L 2 (0) satisfy the properties assumed in the 
previous sections, and A : Q — > M. dxd is symmetric for all x € fi, and uniformly 
positive definite, i.e., there exist constants a, a > such that 

a\£\ 2 <A{x)£,-Z<a\£\ 2 , V x e 0, £ e M d . 

If .4 is piecewise Lipschitz over an initial conforming mesh 7o of SI, i.e., there 
exists CU > such that 

\\A(x) - A(y)\\ 2 < C A \x - y\, Vx,2/eT, VT e To, 

then the convergence results previously presented also hold for this problem. 

Remark 4.11. We have assumed the use of linear finite elements for the dis- 
cretization (see (9)). It is important to notice that the only place where we 
used this is in Proposition 3.3. The rest of the steps of the proof hold regardless 
of the degree of the finite element space. The use of linear finite elements is 
customary in nonlinear problems, because they greatly simplify the analysis. 
The numerical experiments of the next section show a competitive performance 
of the adaptive method for any tested polynomial degree (up to four). 



5 Numerical Experiments 

We conclude this article reporting on the behavior of the adaptive algorithm 
for some particular nonlinear problems. In the first subsection we study the 
convergence rate in terms of degrees of freedom for an exact solution and differ- 
ent functions «(■,■)■ I n the second subsection we show the performance of the 
algorithm when approximating an unknown solution of a prescribed curvature 
equation. 

5.1 Exact solution 

Let us consider the problem 

| -V • [a(\Vu\ 2 )Vu] = / inO 

u = g on <90, 

where C I 2 is the L-shaped domain given in Figure 1. In the following exam- 
ples, in order to study experimentally the behavior of the Adaptive Algorithm, 
we consider (23) with different choices of the function a, defining / and g in each 
case so that the solution of the problem is the function u depicted in Figure 1, 
given in polar coordinates by 

u(r, if) = ri sin ( ^-ip J . (24) 



We consider the Adaptive Algorithm using different marking strategies, 
namely, global refinement, maximum strategy with 9 = 0.7 and Dorfler's strategy 
with 9 = 0.5 (see [11]). 

We implemented the Adaptive Algorithm using the finite clement toolbox 
ALBERTA [11]. We iterated the algorithm until the global error estimator was 
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Figure 1: The domain fi where the problem (23) is posed and the function u 
which is the solution of the problem in each example. 



below 10 -6 or the number of degrees of freedom exceeded 5 x 10 5 . We tested 
the limits of our theory by trying with some functions a which did not satisfy 
all the assumptions of the theoretical results above. 

Example 5.1 (Optimal rate of convergence when a satisfies the hypotheses). As 
a first example, in order to study experimentally the rate of convergence of the 
Adaptive Algorithm, we consider 

which satisfies the hypotheses to guarantee the convergence (see Figure 2), i.e., 
a is a C^-function, and there exist positive constant c a and C a such that 

c a < a(t 2 ) + 2t 2 a'(t 2 ) < C a , Vi > 0, (25) 

and 

a is monotone decreasing, i.e., a'(t) < for all t > 0. (26) 

a(t) a(t 2 ) + 2t 2 a'(t 2 ) 



1.4 
1.2 
1 




°' 5 10 20 30 40 50 °' 2 10 20 30 40 50 



Figure 2: The function a(t) = j^j + \, of Example 5.1 satisfies the properties 
we require to guarantee the convergence. 

In Figure 3 we plot the i? 1 (fi)-error versus the number of degrees of free- 
dom (DOFs), for finite elements of degree £ = 1,2, 3, 4. In this case, the rate of 
the convergence is optimal for adaptive strategies, that is, ||u — m) = 
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0(DOFs fc ). For global refinement, the observed order of convergence is 

— 1/3 

DOFs fc ' for all tested polynomial degrees, due to the fact that the solution u 
belongs to H 1+5 (Cl), for all < S < §, and does not belong to i? 1+2 / 3 (f2). 

Note that, although the theory only guarantees the plain convergence for 
linear elements (cf. Theorem 4.9), the numerical results suggest that the method 
works for any polynomial degree (see Remark 4.11), and the convergence rate 
is optimal. 




Figure 3: Error versus DOFs for Example 5.1. We present the i? 1 (ri)-error 
between the exact solution and discrete solutions, versus the number of degrees 
of freedom (DOFs) used to represent each of them. We note that the conver- 
gence rate is optimal for the considered adaptive strategies, but not for global 
refinement, due to the fact that the solution u is not sufficiently smooth. In 
this case, a(t) = + | satisfies all the properties established to guarantee 
the convergence with linear finite elements. The numerical experiments suggest 
that the method converges with optimal rate for any polynomial degree. 



Example 5.2 (About the hypothesis (25)). We consider the function 

1 1 

a(t) = 1 , t > 0, 

w 1 + f 10 

which is monotone decreasing, i.e., satisfies (26), but not (25), as it is shown in 
Figure 4. Since (25) guarantees the wcll-posedness of problem (23) (uniqueness 
and stability), we could be facing an example with multiple solutions. 

In Figure 5 we plot the H 1 (n)-eivoi versus the number degrees of freedom, 
for different polynomial degrees. For £ = 3 and i = 4 the algorithm stopped 
with an estimator below the desired tolerance 10~ 6 , although the error is around 
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Figure 4: The function a(t) = + of Example 5.2 does not satisfy a(t 2 
2iV(i 2 ) > for all t > 0. 



10~ 2 in all cases. On the other hand, as we can see in Figure 6, the global error 
estimator decreases with optimal rate for the adaptive strategics, indicating that 
the adaptive algorithm may be converging to another solution of the nonlinear 
problem, different from the one given by (24). 

Polynomial degree i = 1 Polynomial degree 1 = 2 




Optii 



um Strategy 
Strategy 



D O F ■ 



Optii 



um Strategy 
■ Strategy 




Figure 5: Error versus DOFs for Example 5.2. We present the 7? 1 (ri)-error 
between the exact and discrete solutions, versus DOFs. We observe that the 
method does not converge, but the error stagnates around 10~ 2 . Since a(t) = 
+ does not satisfy the condition which guarantees uniqueness of solutions, 
and based on the fact that the a posteriori error estimator does tend to zero 
(see Figure 6) we conclude that the method converges to another solution of 
the same problem, different from the one given by (24). 
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Based on this remark, it seems that the adaptive algorithm converges to a 
solution u\ such that ||it — UiHifim) rj 1CP 2 . We recall that a does not satisfy 
the condition (25) which guarantees the uniqueness of solutions. 





Figure 6: Estimator versus DOFs for Example 5.2. We present the global 
error estimator rjk(Tk) versus DOFs used to represent each discrete solution. 
We note that for the adaptive strategies the global error estimator decreases 
with the optimal rate for the H 1 -error, although the error does not tend to zero 
(see Figure 5). In this case, a(t) = + ^ does not satisfy the condition 
which guarantees uniqueness of solutions. It seems that the method converges 
to another solution of the problem. 



Example 5.3 (About the hypothesis (26)). We consider 

1 3 

a(t) =--e"5*+l, t>0, 

which satisfies the hypothesis (25) related to the well-posedness of the problem 
but not (26), because a is monotone increasing (see Figure 7). 

In Figure 8 we plot H 1 (fl)-error versus the number of degrees of freedom, 
for finite elements of degrees I = 1,2,3,4. Note that in this case the optimal 
convergence rate DOFs - ^ 2 is still observed for the adaptive algorithm. This is 
an indication that the assumption (26) about a being monotone decreasing can 
be superfluous, and only an artificial requirement for the presented proof (see 
Lemma 4.1). A more detailed study about this hypothesis is subject of future 
research, and beyond the scope of this article. 
Example 5.4. Finally, we consider an extreme case, with 

a(t)=2 — T , t>0, 
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Figure 7: The function a(t) = 
esis (25) but is monotone increasing 



i -±t 

- e 2 r - 



-1 of Example 5.3 satisfies the hypoth- 
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Figure 8: Error versus DOFs for Example 5.3. We present the H 1 ^) -error 
between the exact and discrete solutions, versus DOFs. We note that the con- 
vergence rate is optimal for the considered adaptive strategies, although the 
function a(t) = — ^ c~ a* +1 is not monotone decreasing. This could mean that 
this hypothesis is not necessary for the convergence of the Adaptive Algorithm, 
which performs better than expected by our theory. 
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which satisfies (25) and (26), but lim t _ >0 + a '{t) = — oo, as can be observed 
in Figure 9. This means that a is not Lipschitz continuous. Since we only 
require that |ta'(i)| is bounded (cf. (25)), it still satisfies the assumptions of 
the convergence theory, and optimality is observed regardless of the fact that 
lim t _j. + a '(t) = —oo- 



c(t) c(t 2 ) + 2t 2 a'{t 2 ) 




Figure 9: The function a(t) = 2 — y+v? °^ Example 5.4 satisfies our assumptions 
which do not require that a is Lipschitz continuous. 



Figure 10 shows if 1 (ri)-error versus DOFs, for polynomial degrees I = 
1, 2, 3, 4. We obtain optimal convergence rate DOFs~ £ / 2 for the adaptive strate- 
gies. Thus, even when a is not Lipschitz, the convergence rate is optimal. As 
a consequence we conclude that it should not be necessary to make additional 
assumptions in order to prove optimality. 



5.2 Unknown solution 

In this section we use the Adaptive Algorithm to approximate a solution to a 
prescribed mean curvature problem. We consider the problem 



V 1 + iv«| s 



with 



/(*) 







= / in = (-1,1) x (-1,1) 
u = on <9fl, 

if \x\ < 1/3, 
3 if 1/3 < |x| < 2/3, 
otherwise. 



(27) 



The function a(t) = 1/ VI + t corresponding to this equation does not fulfill (25) 
because a(t),a'(t) — > as t — > oo. Nevertheless, for many right-hand side func- 
tions /, as the one stated above, the solution belongs to W 1,ao {£l). This implies 
that |Vu| is bounded and a could be replaced by a function satisfying (25) 
without changing the solution. This is not needed in practice, but is rather a 
theoretical tool for proving that the assumptions hold in some sense. 

We experimented with several right-hand sides / and observed that the 
method performs robustly whenever u € W ,oo (0). When the solution has an 
unbounded gradient, the method produces a sequence with a maximum value 
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Figure 10: Error versus DOFs for Example 5.4. We present the 7? 1 (^l)-error 
between the exact and discrete solutions, versus DOFs. We note that the con- 
vergence rate is optimal for the considered adaptive strategies, even though the 
function a(t) = 2— has an infinite derivative at 0. This example falls inside 
the theory, and the fact that a is not Lipschitz continuous does not destroy the 
optimality of the sequence of discrete solutions. 
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increasing to infinity. This is a drawback of our method, since it cannot be 
used to approximate singular solutions (not belonging to W 1,00 (Q)), if a does 
not satisfy (25). It is worth recalling that adaptivity is still a good tool for ap- 
proximating regular solutions, specially taking into account the fact that in this 
context regularity is a concept relative to the polynomial degree. Using adap- 
tivity with higher order finite elements can improve drastically the performance 
even when the solution is in H 2 (Q) n W 1,0O (£i). 

We present in Figure 11a picture of the solution obtained with our method 
and several meshes at different iteration steps, for polynomial degrees 1, 2 and 
3. It is worth observing the stronger grading obtained for higher polynomial 
degrees. This obeys the fact that we mention above, that the singularity of the 
solution is relative to the polynomial degree of the finite element space. Recall 
that in practice adaptive methods with polynomials of degree i on domains in 
R d converge with order DOF~ £ / d . In order to obtain this rate with uniform 
meshes, we would need the solution to belong to H e+1 (£l). 
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